!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!sav2008 pn=fn(r)
pn=fn1(r,fnr)
pn=fn2(r,fnr,fnrr) !sav2008
!!!!!!!read data !!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!sign of driven current in right coordinate system {dro,dteta,dfi}:
!!!curdir=+1.0 for current drive in positive direction "dfi"
!!!curdir=-1.0 for current drive in negative direction "dfi"
!!!spectrum Nz>0 is along dfi>0 and Nz<0 is along dfi<0
!!!it is also OK if Npar is used instead of Nz, but for Btor>0, that is along dfi>0
curdir=-dble(ispectr)
!!!!!!!!!!!!! begin iterations !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
find nevyazka
!!!!!!! prepare to the next iteration !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! if(q_cond.le.pabs0.and.pchg.lt.pgiter)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!! and multiplied by 10^-10 in dfind() !old dconst=vt0/(cnyfocpmepn2) ! dj(i)=dql(i,j)dconstvto !D_normir ! call linf(vrj,dj,vij(i,j),ddout,klo,khi) if(ddout.le.1.d0) ddout=zero
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(Spectrum) | :: | spectr | ||||
| real(kind=8), | dimension(*) | :: | outpe | |||
| real(kind=8) | :: | pe_out |